I©  RDECOM 


TECHNOLOGY  DRIVEN.  WARFIGHTER  FOCUSED. 


Math-Based  Simulation  Tools  and  Methods 

Sudhakar  Arepally,  CRSR  Team,  US  Army  RDECOM-TARDEC 
10  October  2007 


IM 


HMMWV  30-mph  Rollover  Test 
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Soldier  Gear  Effects 
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Occupant  Performance  in  Blast  Effects 


Joint  Constraint  Force-1 
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Anthropomorphic  Test  Device 
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Figure  2.1  Hvbjsd  EH  5C*  percentile  dummy  modej:  elUptc-id  mode',  (left;  and  facet  model 
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Lower  Limb  Model 


Brain  Model 


Abdomen  Model 


Thorax  Model 


Wayne  State  -  ANSIR 


Models 


MADYMO-  FACET 
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MADYMO  FE-  HUMOS2 


a.)  soft  tissues  b)  sketetrn  t)  orsanj 
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Rigid  Body  Modeling 


KINEMATICS  OF  A  RIGID  BODY 


Position  vector  of  P,  a  point  on  body  i  in  reference  space  is  given  by, 


A  multibody  system  is 
a  set  of  bodies 
interconnected  by  a 
kinematic  joints 


Xl  =rl+COixxi,  where  (Oi  is  the  angular  velocity  vector  of  body  i 
And  2nd  time  derivative  equals, 

X.  =  r.  +  d).  x  x  +  (0-  x  (co,  x  X.. ) 

L  L  L  L  L  'I  v  ' 
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Rigid  body  modeling  -  continued 


KINEMATICS  OF  BODIES  CONNECTED  BY  JOINT 

—j~  —i  —ij  —ij  —  Ji  ^  1 

Where, 

Cg  =  time  dependent  direction  cosine  matrix,  joint  coordinate  system  on  body  i 


relative  to  body  local  coordinate  system  on  body  i 

D  -  =  direction  cosine  matrix,  orientation  of  joint  coordinate  system  on  body  j  with 
relative  to  joint  coordinate  system  on  body  i 

A-  =  direction  cosine  matrix,  orientation  of  local  coordinate  system  of  body  i 
A  .  =  direction  cosine  matrix,  orientation  of  local  coordinate  system  of  body  j 


rj=ri+cij+dij~cji  =-2 

fj  =  position  vector  of  the  origin  of  the  local  coordinate  system  on  body  j 


r. •  =  position  vector  of  the  origin  of  the  local  coordinate  system  on  body  i 

c  =  position  vector  of  the  origin  of  the  joint  coordinate  system  on  body  i  relative  to  the  origin 
V 

of  the  local  coordinate  system  of  the  corresponding  body 

c  =  position  vector  of  the  origin  of  the  joint  coordinate  system  on  the  body  j  relative  to  the  origin 

fi 

of  the  local  coordinate  system  of  the  corresponding  body 

dH  =  vector  from  the  origin  of  the  joint  coordinate  system  on  body  i  to  the  joint 

coordinate  system  on  body  j 

Taking  first  time  derivatives,  yield  the  following  expressions  for  angular  and 
linear  velocities  - 


0)  .  =  CD.  +0)..  ^3 

J  I  u 


Two  bodies  interconnected  by  an 
arbitrary  kinematic  joint;  motion  of 
body  j  is  described  relative  to  parent 
body  i 


rj=ri+a)iXCij+dij-(OjXCji 


=>  4 
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Tnrtsomc 


Equations  of  motion  (Newton  -  Euler)  of  rigid  body  i  referred  to  its 
center  of  gravity  are : 


J:  6);  +6);  X  J,  O);  =  T: 


Where, 
mj  =  mass 

Jt  =  inertia  tensor  w/r  to  center  of  gravity 

(Ox  =  angular  velocity  vector 

Ft  =  resultant  force  vector  relative  to  the  CG 

Tt  =  resultant  torque  vector  relative  to  the  CG 

Ft  and  Tt  both  include  constraint  forces  and  torques  due  to  joints 

Using  principle  of  virtual  work  and  variation  of  position  vector,  8ri  and 

a  variation  of  orientation,  8k l ,  the  resulting  equations  summed  for  all 

bodies  in  the  system : 
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Rigid  body  modeling  -  continued 


Expression  for  2nd  derivative  of  the  joint  degrees  of  freedom  are  obtained : 


<?.,  =KijL+Q,.: 


—V 


— lJ 


Where, 

Y .  is  a  6  x  1  column  matrix  that  contains  components  of  linear  and  angular 

acceleration  of  the  coordinate  system  of  the  parent  body  i 

n..  x  6  matrix  M..  and  the  n..  x  1  column  matrix  of  Q  depend  on  the  inertia  of 

j  j  j  - ij 

the  bodies  and  the  instantaneous  geometry  of  the  system.  Q  depends  additionally 

— U 


on  the  instantaneous  geometry  of  the  system  and  the  applied  loads. 
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Rigid  body  modeling  -  continued 


Tnrtsomc 


Numerical  integration  methods : 

The  equations  of  motion  form  a  system  of  coupled  non  -  linear  second 
order  differential  equations : 


q  =  h\q,  q,  t  J,  with  initial  values  of  q0  and  q0 


q  =  column  matrix  with  generalized  coordinates;  joint  position  degrees  of  freedom 
q  =  joint  velocity  degrees  of  freedom 
q  =  joint  acceleration  degrees  of  freedom 

Equations  of  motion  are  solved  numerically  using  the  available  three  methods : 

-  a  modified  Euler  method  with  a  fixed  or  variable  time  step 

-  a  Runge  -  Kutta  method  with  a  fixed  or  variable  time  step 

-  a  Runge  -  Kutta  Merson  method  with  a  variable  time  step 
Modified  Euler  method : 

One  step  method  with  fixed  time  step  ts;  second  order  differential  equations  are  integrated 
explicit  Euler  method  which  gives  velocity  variables  at  time  point  tn+l  =tn+  ts 


The  velocity  variables  are  integrated  using  the  implicit  Euler  method  which  gives 
solution  for  the  position  variables  at  time  point  tn+l , 
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Finite  Element  Methods 


TnnsoiEi 


Displacement  Matrix 

l<5f  —  \u.VM  ,v  m  , v  r 

L  Je  l  i 5  i ?  j  ?  j ?  m  ’  mJ 

or,  displacement  function  at  any  point  within  the  element  is 
Stress,  Strain  and  Elasticity  Matrices 

W=M4 

Formulation  of  Finite  Element  Methodology  using  the 
principle  (minimization)  of  total  potential  energy  of  the  system  - 

tlfawi-towity-iitonn*-* 

i  i  4 

or,£{A<JM4{4-{eU 

1 

the  stiffness  matrix,  \k\ ,  and  nodal  force  matrix,  \Q\ , 

(due  to  body  force,  initial  strain,  and  surface  traction)  are 

fol  =  l  ptf  {FW  + 1  [fif  Me.. }dV  +  j>f  {T}ds 

And  the  system  of  equations  for  equilibrium  for  nodal  forces  for  the  entire  structure, 
where, 

W  =  E[*Land{e}=E{e}« 


>  y,  v 


Typical  finite  element  ‘e’ 
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Injury  Criteria  development 


Scalp,  Skull,  Meninges  and  Brain 


SPACE 
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Injury  Criteria  development  -  continued 


Head  Injury  Criteria : 

1  S 

HIC  =  max 

1  2 

- [  a(t)dt 

^2  —  h  tx 

A*  m  w/ 

(c  — 0 

Where,  tx  and  t2  are  two  arbitrary 
times  during  acceleration  pulse 

HIC  is  limited  to  36  ms 

Lognormal  distribution  (pdf) : 


/(*”)= 


1 


<7,  a/2/t 

T 


e 


t-h 

a  . 

V  T  J 


The  probability  of  skull  fracture  (AIS^2)  is  given  by  tjlie  formula 

(  h\(  HIC)  -  /  A 


p  ( fracture )  =  N 


(7 


TRf7D£=l 

Wayne  State  Tolerance  Curve 


where  N( )  is  the  cumulative  nonnal  distribution.,  p  =  6.96352  and  <5=  0.84664. 


Injury  Risk  Curve  for  HIC 
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Non  -  linear  Constrained  Optimization  Problem : 
Minimize :  F  (v)  Obj  ective  function 

Subject  to : 


g,(x)<0  j  =  \,m 

hk{x)  =  0  k  =  \,l 


X'<XS  <X“  i  =  l 


n 


X , 


X2 

^3 


where,  X  =  < . 


Inequality  constraint 
Equality  constraint 
Side  constraints 


ivj 

Necessary  and  sufficient  conditions  for  global  optimum 
for  unconstrained  problems  - 
VF(x)  =  0  (the  gradient  must  vanish) 


itimization 


TnnotEC 
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Optimization  -  continued 


where, 


d 


F{X) 


d 

dX2 


F{X) 


VF(X)  = 


d 


F(X) 


and  the  Hessian  matrix  must  be  positive  definite 


H  = 


32F(x)  32F(x) 

ax,2  ax,ax2 

d2F  a2F(x) 
ax2ax,  ax2 


a2F(x) 

ax,ax2 

a2F(x) 

ax^x; 


a2F(x)  a2F(x)  a2F(x) 
ax„ax,  ax„ax2 ax2 


And  for  constrained  problems,  if  vector  X*  defines  the  optimum  design, 
Kuhn  -  Tucker  conditions  are  necessary  and  sufficient : 

1 .  X*  is  feasible 

2-Ajgj(x")=0  j  =  \,m  Aj>  0 

3  .VFfx* )+  £  A.Vg.  fx* )+  X  4+tVAt  (x* ) 

7=1  &=1 

X>0 
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•  Crashworthiness  and  occupant  protection  in  motor  vehicles 
comprise  the  following  significant  mathematics  applications  - 

-  Matrix  operations 

-  Ordinary  and  partial  differential  system  of  equations 

-  Lagrangian  operations 

-  Fourier  transforms 

-  Taylor  series 

-  Finite  difference  methods 

-  Implicit  and  explicit  finite  element  methods 

-  Statistical  methods  -  probabilistic  and  regression  analysis 
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